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Abstract 

We construct the family of spin chain Hamiltonians, which have affine 
quantum group symmetry U q g. Their eigenvalues coincide with the eigen- 
values of the usual spin chain Hamiltonians, but have the degeneracy of 
levels, corresponding to affine U q g. The space of states of these spin chains 
is formed by the tensor product of fully reducible representations of quan- 
tum group. 

The fermionic representations of constructed spin chain Hamiltonians 
show that we have new extensions of Hubbard Hamiltonians. All of them 
are integrable and have affine quantum group symmetry. The exact ground 
state of a such type model is presented, exhibiting superconducting behav- 
ior via 77-pairing mechanism. 
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1 Introduction 



The interest to strongly interacting electron systems was raised after discovery 
of high-T c superconductivity [|T| , where it is believed that the interplay between 
magnetism and insulating behavior has a crucial role. The Hubbard model [0 
is the simplest example of strongly correlated electrons and was solved in one 
dimension via Bethe ansatz technique ||. Subsequently numerous publications 
appeared investigating one-band (multi-band) Hubbard models in one and two 
dimensions, most of which used approximate or numerical methods and only few 
exact results are known. One of these exact results based on so called //-pairing 



mechanism of superconductivity, the idea of which was introduced by Yang WA 



15|1 for the Hubbard model. It appeared, that there is certain type of states 
exhibiting off-diagonal long-range order (ODLRO), the concept, which had been 
developed in 0]. In 0, |5], || was shown that ODLRO implies Meissner effect and 
flux quantization and thus can be regarded as definition of superconductivity. But 
ground state of simple Hubbard model is not superconducting j| . The first model 
with truly superconducting ground state of 77-pairing type is the supersymmetric 
Hubbard model, introduced in [|I6|, |T7||. This model is a certain extension of 



Hubbard model by additional nearest-neighbor interactions and some restriction 
on hopping mechanism. 

It appeared that ^-pairing mechanism is not an exotic phenomenon and re- 
cently some extended Hubbard models were defined, superconducting ground 
state of which could be found exactly [0, [8], || [II], [13| . The main property of 
these extensions is the fact, that they besides the ordinary electron hopping and 
Hubbard interaction terms contain also bond-charge interaction, pair-hopping, 
nearest-neighbor Coulomb interaction and XXX spin interaction terms. 

In this article we develop the new technique and define a family of other 
extensions of one dimensional Hubbard Hamiltonians, one of which exhibits su- 
perconductivity based on ^-pairing. The defining property of this family is the 
fact, that they have U q sl{2) affine quantum group symmetry. 

Quantum group symmetry plays essential role in integrable statistical models 



21| , |22|| and conformal field theory [^, [24], [25] 



It is well known that many integrable Hamiltonians have quantum group 
symmetry. For example, XX Z Heisenberg Hamiltonian with particular boundary 
terms is [7^ ^-invariant p4 |. The infinite XX Z spin chain has larger symmetry: 
affine U q sl2 pB| . The single spin site of most considered Hamiltonians forms an 
irreducible representation of Lie algebra or its quantum deformation. 

In (IIJ we have constructed the family of spin chain Hamiltonians 7i, which 
have affine quantum group symmetry U q g, using intertwining operators between 
tensor products of its spectral parameter dependent reducible representations. 
These intertwining operators commute with U q g by definition, whereas Hamilto- 
nians 7i commute with U q g by construction. The space of states of these spin 
chains is formed by the tensor product of the fully reducible representations. We 
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have shown that the model, considered in P7|, which corresponds to some ex- 



tension of the Hubbard Hamiltonian in the strong repulsion limit, is a particular 
case of our general construction. The affine quantum group symmetry leads to a 
high degeneracy of energy levels. 

The energy levels of these spin chains are formed on the states, constructed 
from the highest weight vectors of quantum group representations. In particular 
cases the restriction of the considered spin chains on these states gives rise to 
Heisenberg model or Haldane-Shastry long-range interaction spin chain. Hence, 
we have the generalizations of these Hamiltonians, which have affine quantum 
group symmetry. 

It appeared that the fermionization of the simplest examples of these spin 
chain Hamiltonians gives the family of extensions of Hubbard model. 

In the sec. 2 we give the definitions of quantum Kac-Moody group U q g and 
properties of i?-matrix. In sec. 3 we develop the technique for construction of 
nearest-neighbor spin chain Hamiltonians, which have U q sl2 quantum group sym- 
metry. In sec. 4 we perform this construction to obtain Haldane-Shastry type 
long-range spin chain Hamiltonians. 

In sec. 5 we consider the fermionization of previously constructed integrable 
spin chains for some simplest particular cases, in particular, when the site is 
four dimensional and there is no dependence on quantum group deformation 
parameter q. We obtain in such way some integrable extensions of Hubbard 
model. In the case, when the space of states of each site is a direct sum of spin-1 
and trivial multiplets, one obtain Hubbard model with doping, pair hopping and 
nearest-neighbor interaction of holes. 

In the case, when the space of states of each site is a direct sum of two 
two dimensional representations, one obtain Hubbard model with bond-charge 
interaction, density-density interaction and 'boson-boson' interaction between 
nearest-neighbor hole and double occupied site. We find the ground states in the 
fixed particle number sectors and phase space region, where this model exhibits 
//-pairing superconductivity. In the particular case, when the amplitude of last 
two interactions vanishes, one obtain well known extended Hubbard model with 



bond-charge interaction O 



2 Definitions 



Let us recall the definition of quantum Kac-Moody group U q g (U q g) |20], ^T], ^2| . 
It is generated by elements e,, hi satisfying the relations 

[/ij,ej] = Cijtj [hi, fj] = —Cijfj (2-1) 
[eijj] = 5 lj [h} q 

and g-deformed Serre relations, which we don't write. Here i — 0, . . . , n for U q g 
and i — 1, . . . , n for U q g, q is a deformation parameter, [x] := {x q — x~ q )/(q— q^ 1 ), 
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Cjj is a Cartan matrix of corresponding affine Lie algebra g (finite Lie algebra g) . 
On U q g (U q g) there is a Hopf algebra structure: 

A(e;) = ki ® ei + e t ® fc^ 1 A(fcf x ) = fcf 1 ® fcf 1 

A(/<) =h®fi + fi® K 1 (2.2) 

where fcj := q~ . This comultiplication can be extended to L-fold tensor product 
by 

L 

A L_1 (ei) = ® . . . ® h ® OA;^ 1 ® . . . ® A;," 1 

L 

A L_1 (/i) = ^ fei (8 . . . ® A* <g> /» &K 1 ® . . . ® A^ 1 

A^ifef 1 ) = kf 1 <g> . . . ® Jfef 1 
There is an opposite comultiplication A, which is obtained from (|2.2|) by replacing 

For general g the representations of quantum group [/ g <7 are is one to one 
correspondence to the representations of nondeformed Lie algebra g. Denote by 
Va the irreducible [/gg-multiplet with highest weight A. It possesses a highest 
weight vector v x , such that 

e iV ° x = tnv° x = \(hi)v x , i = l,...,n (2.3) 

Let g be an affine algebra and g be its underlying finite algebra. Then for any 
complex x there is the g-deformation of loop homomorphism p x : U q g — » U q g, 
which is given by [^0] 

Px(e ) = xf e p x {fo) = x~ l e e p x (h ) = -he 
px(ei) = ei Px(fi) = fi Px(hi) = hi, 

where i = 1 . . . n and 9 is a maximal root of U q g. Using p x one can construct 
the spectral parameter dependent representation of U q g from the representation 
of U q g. 

Let ~Vx( x i) an d ^2(^2) are constructed in such way irreducible finite dimen- 
sional representations of U q g with parameters x\ and x 2 correspondingly. The 
U q ^-representations on Vi(xi) <g> V 2 (x 2 ) constructed by means of both A and A 
are irreducible, in general, and equivalent: 

R(x 1 ,x 2 )A(g) = A(g)R(xi, x 2 ), Vg E U q g 

The i?-matrix R{x-i : x 2 ) depends only on X\jx 2 and is a Boltsmann weight of 
some integrable statistical mechanic system. 
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3 Quantum group invariant Hamiltonians for re- 
ducible representations 

Let V = (B^LiV^ is a direct sum of finite dimensional irreducible representations 
of U q g. We denote by V(xi, . . . , x^) corresponding affme U q g representation with 
spectral parameters xf. 

N 

V(x 1 ,...,x N )=®V Xi (x i ) (3.1) 

i=l 

We consider the intertwining operator 
H(x 1 , ...,x N ): 

V(xi, . . . , x N ) <g> V(xi, . . . , x N ) V(xi, . . . , x N ) ® V(xi, . . . , x N ), 
[H(x 1 ,...,x N ),A(a)) = 0, VaeU q g 

If V = V\ consists of one irreducible component then if is a multiple of identity, 
because the tensor product is irreducible in this case. To carry out the general 
case let us gather all equivalent multiplets together: [] 

M 

V(x 1) ...,i JV )=0N A ,®V A ,(i ! ), (3.2) 

i=l 

where all Va( M nonequivalent irreps and ~ C^* have a dimension 

equal to the multiplicity of V\ t (xi) in V(xi, . . . ,xn)- Note that J2 i=1 = N. 
By the hat over the tensor product we mean that U q g doesn't act on Na^Va^x^) 
by means of A but acts as id® g. 
So, we have: 

V(xi, . . . , x M ) ® V(xi, . . . , x M ) 

= (®fiiN Al ® V Al (x,))(g)(©^ 1 NA I ®VA l (x i )) 

M 

= Na^Na^^a^xO^Va^x,)) 

Now, V Ai (xj) ® V Aj (xj) is equivalent only to itself and to V Aj (xj) (g) V\X x i) (f° r 
i j) by applying the intertwining operator R(xi/xj) = PR(xi/xj), where P is 
tensor product permutation: P(vi®v 2 ) = v 2 ®Vi. So, the operator H(xi, . . . , %), 
commuting with U q g on V(x\, . . . , xm) <8> V(xi, . . . , xm) has the following form: 

H\-N Xi ®H\ j ®Vx i ®Vx j = ^i® id v Ai ®v Aj + B ij ®>Rv x .^Vx j ( x i/ x j) ( 3 -3) 

1 The [/^-equivalence of V,^ (xj) requires that the spectral parameters Xi and highest weights 
Xi are the same. 
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where Ay and Bij are any operators on Na^©Na.. Note that in one special case 
the formula (|3.3| ) simplifies. Consider the situation if there is only one nontrivial 
representation Va, i.e. A 7^ 0, in the decomposition ( |3.2| ) of Y(xi, . . . ,xn). As 
,g,Vo = -Rv ®v A = id, the intertwining operator ( jOD doesn't contain spectral 
parameters Xi and deformation parameter q. So, H is invariant under the action 
of quantum group for all values of deformation parameter. Of course, this is much 
larger symmetry and it can be written in another form without using of q. 

To write down the action ( j3.3|) more explicitly we introduce the projection 
operators 

XS=\a){b\, (3.4) 

where the vectors \a) span the space V. In accordance with the decomposi- 
tion ( |3.2j ) we use the double index a = (n», Oj), i = 1,...,M where the first 
index n« = 1, . . . , Ni characterizes the multiplicity of Va, and the second one 
a>i = 1, . . . , dimVA^ is the vector index of Va,. Then the formula ( |3.3| ) can be 
represented as 

H(A B)=jt[ £ A^ £ ® X^;% (3.5) 

di,dj 

+ E £ Ri^i/x^X^®X^% 

where B u = and we use the matrix form of .R-operator: 

#v Aj ®v* (>i/%)k) <g> \aj) = R ijTJ.( x i/ x i)\ a 'i) ® K>- 

a. ,a. 

Let us consider some particular cases of this general construction. 

(i) Let V(x) = V(x,x) = V\(:e) © ~V\(x). The second term in (|3.5|) is absent 
in this case and H has factorized form: 

H = A®id Vx ® Vx , A = A%™ (3.6) 

where ni, 712, rn\, m-i = 1,2 are indexes, which mention each Va- In the matrix 
form ( |3.6|) can be written as 

h= E ^E4^®4^ ( 3 - 7 ) 

n,i,rij ,mi,rrij o,i,a,j 

(ii) Let now V(xi,x 2 ) = V Al (xi) © Va 2 (x 2 ) (V\ t (xi) are mutually nonequiv- 
alent). Then H acquires the following block-diagonal form 



H(x 1 ,x 2 ) 



( a ■ id \ 

c • id d ■ ^12(^1/^2) 

e • ^21(^2/^1) /'id 

V 0-id J 



(3i 
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Here to be short we used the notation 

R 12 (xi/x 2 ) = R Vxi ^ Vx2 (xi/x 2 ), R 2 i(x 2 /x 1 ) = R Vx2 ® Vxi (x 2 /x 1 ). 

(hi) Let us choose g = sl(2) and V = V S ©V , where V s is 2s + 1-dimensional 
spin-s representation of U q sl 2 and Vq is the trivial one dimensional representation. 



This case was considered in 27 



Following ]27j from the operator H the following Hamiltonian acting on W = 
\® L can be constructed: R 

u L-l 

H = £ Hii+i (3.9) 
i=i 

Here and in the following for the operator X = J2i x i ® yi on V <g> V we denote 
by its action on W defined by 

Xij = id ® • • • ® id <g> <g>id . . . <g> id <g> ^ <g>id <g> . . . <g> id (3.10) 

' * j 

By the construction, 7i is quantum group invariant: 

[H,A L - 1 (g)] = 0, WgeU q g 

Let us define the projection operator V\ on V for each class of equivalent 
irreps (A i; x^), i — 1, . . . , M: 

$>i = id, ^ 2 = ^ 

i=l 

Their action on the tensor product W is given by 

L 

Pi = id <g) . . . id ® "Pi <g>id . . . <g> id 

It is easy to see that these projections commute with Hamiltonian 7i and quantum 
group U q g: 

[V u H} = 0, [Pi,U q g] = (3.11) 

Denote by W pi „ _ PM the subspace of W with values Pi of V{ on it. Then we have 
the decomposition 

W W P1 ... PM (3.12) 

Pl.---.PA/ 
pi + ...+PM=L 



J Here and in the following we omit the dependence on xi 
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Let V° is the linear space, spanned by the highest weight vectors of U q g- 
module V: 

V° := 

where v®. G V A . is highest weight vector ( f2 . 3| ) . We define also W° := V° ® L . 
The space W° is 7i- invariant. This follows from ( p.3|) . For general q the action 
of U q g on W° generate whole space W. Indeed, the [/^-action on each state of 
type Vy <8> • • • <8> v A . generates the whole space V A . <S> ■ ■ ■ <8> V A . , because the 
tensor product of finite dimensional irreducible representations of affine quantum 
group is irreducible |28[ . 

Consider now the subspace W° = W° n W P1 ... PM . We have the decom- 
position, which is inherited from ( |3.12| ) 

w °= w L. P m ( 3 - 13 ) 

P1.—.PM 
Pl+...+PAf=L 

Note that 

d p ,.. PM := dimWj^ = ( ^ L ^ ) iVf . . . iVlr 

Let us define by ifo an d Ti-o the restrictions of H and 7f on V° <S> V° and 
W° correspondingly: H := if| v o<g)V ; := 7~L\w°- It follows from ( |3.11| ) 



that Hamiltonians 7i and 7io have block diagonal form with respect to the de- 
compositions ( 3.12 ) and (|3.13| ) correspondingly. Every Hamiltonian eigenvector 



w° G W" „ with energy value E ari „ gives rise to an irreducible U a q- 
multiplet W api p of dimension 

Af 

dimW api ... PM = n(dimV Afc r (3-14) 

fc=i 

On W„ the Hamiltonian H. is diagonal with eigenvalue E n . In partic- 
ular case when all V Ai are equivalent, the degeneracy levels are the same for all 



E ap and are equal to (dimV A ) L . Note that 



a P1 ... PM =l 

P! + ...+p M =L 



m—PM 

dimW dimW a P1 . ,, u 



Af M 

Pk 



pi + ...+PM=L 

(j> fc dimV Afc J 
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according to the decomposition 



fl.--.PM a P1 ... PM = l 



•"P1---PM 

y„-, r,,,=l 

pi + ...+p M =L 

Now, we suggest that one can obtain the energy spectrum E api _ PM for 7Y - 
Then for the statistical sum we have 

Zn ((3)= E E exp(/3£ aj>i ... PM ) (3.15) 



P1.--..PM a Pl ... PM =l 
pi+...+p A /=L 



Then the statistical sum of 7i has the following form: 

M d vl ... PAI 

Zn(P)= E II(dimV Afc r E exp(/3£ ajJi ... p J (3.16) 

Plv.PAf fc=l a P1 ...p A/ =l 
Pi+...+p A/ =L 

So, if the underlying Hamiltonian Tio is integrable and its eigenvectors and eigen- 
values can be obtained, then this problem is solved for 7i also. Performing the 
quantum group on all eigenvectors of an energy level of Tio one can obtain the 
whole eigenspace TL for this level. 

Consider now some particular cases of our general construction. If we choose 
two equivalent representations (the first example above), then dimV = 2 and 
there is one term in decomposition (|3.13|) . H now is the most general action 
on V° g) V°. As a particular case, the XYZ Hamiltonian in the magnetic field 
can be obtained. The correspondence between the spectrum of H, and Tio is the 
simplest in this case because the additional degeneracy of all energy levels is the 
same and for the statistical sums (13.151), (13.161) we have 



ZniP) = (dimV A ) L Z Wo (/3) = (dimV x ) L Z xxz ((3) 
Consider now the second example of i?-operator above. Let us choose for 



parameters in (|3] 



a = g = e = d = l c = f = (3.17) 
One can renormalize i?-matrices in ( |3.8| ) to satisfy the unitarity condition 

Ri 2 (z)R 2 i{z- 1 ) = id 
. Together with ( |3.17| ) this leads to 

H( Xl ,x 2 ) 2 = id (3.18) 
The restriction of Ti. on W° coincides with the XXX Heisenberg spin chain 

^0 = H X XX = E P H+1 = ^ E ( X + <^+l) 

i i 

The space W° ip2 , corresponds to all states with s z = p\/2 value of spin projection 
S z = 1/2^1^- If we return to TC, the energy level degeneracy of each eigenstate 
with the spin projection s z is multiplied by (dimVA 1 ) 2Sz (dimVA 



\L-2s z 
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4 Generalization to long-range interaction spin 
chains 



Let us consider now the generalization of the construction above for the case of 
long-range interacting Hamiltonians. 



Recall the Haldane-Shastry spin chain is given by 29, 30, 31 



< 5 = E^-^> (4-1) 

Here the spins take values in the fundamental representation of sl n and per- 
mutes the spins at z-th and j-th positions. It is well known that the Hamiltonian 
(|4. 1|) is integrable if di has one of the following values 

{j, rational case 

(1/a) sinh(aj), a & R, hyperbolic case 
(L/n) sm(nj /L), trigonometric case 

The trigonometric model is defined on periodic chain and the sum in (|4.1|) is 
performed over 1 < i,j < L. Rational and hyperbolic models are defined on 
infinite chain. 

For simplicity we consider here sl 2 Haldane-Shastry model only. This means 
that the space V° should be two dimensional, so we should have two terms in 
the decomposition ( |3.1| ). 

Now our goal is to construct a long-range interaction model, which on the 
highest weight space W° coincides with ( |4.1|) . So, we take V = V Al © V Aa . This 
is the second example considered in sec. 3. To achieve an permutation action 
on the subspace V° <g) V° let us impose the conditions (|3.17|) on parameters of 
if -operator ( |3.8| ). One can try to generalize the Hamiltonian ( [4.1| ) in this way by 

n-/HS _ \ - 1 tt 

i<j a i~j 

But it is easy to see that then TC HS isn't invariant with respect to quantum group 
U q g. This is because the equation 

R ij (x u x 2 )A L - 1 (g) = A L - 1 (g)R ij (x 1 ,x 2 ), g G U q g 

is valid only for i — j ± 1 . 

To overcome this difficulty let us use instead of the operator Q 

= G[ij\Hj-.\jGy i -Y where = H ii+1 H i+li+2 . . . Hj_ 2 j-i (4.2) 



3 To be short, we omit Xi, a; j -dependence of operators H, G, F-operators in this section. 
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We remark that Fm and G^] act nontrivially on all indexes i, i + 1, . . . ,j. So 
we include them into bracket do not confuse with the definition (|3.10|) . The 
'nonlocal' term like appeared earlier in the construction of quantum group 
invariant periodic spin chains as a boundary term [|3!| [35| . 

The operators Hu + i satisfy the permutation group relations ( |3.18| ) and 

Hi-uHu+iHi-u = Hn + iHi-iiHu + i (4-3) 

These relations obey on W°, because the restriction of Hij on it gives rise to 
permutations P^- 



Hi 



p. 



' ij I Wo ± ij 

On the whole space W the equations ( |4.3| ) are continued using the affine quantum 
group symmetry of H iii+1 . In contrast to the standard realization by P^, the 
relation 

P'i—liPii+lPi—li Pi—li+1 

(4.4) 

isn't fulfilled. It is easy to see from ( |4.2|) and (|4.4| ) that 



[ij]\W° 



P. 



So, the spin chain defined by 



n 



HS 



E 



" d 2 ■ 

i<j 



' [ij] ■ 



(4.5) 



is quantum group invariant and its restriction on the space W° it coincides with 
the Haldane-Shastry spin chain (|4.1| ). 

In some cases nonlocal expression for ( f4.2| ) becomes depending on sites 



V , V 



A 2 



V A ), then 



H 



(4.6) 



i and j only. For example, if is trivial (V^ 1 
-^12(^1/^2) = -^21(^2/^1) = id and if-matrix Q3.8| ) has the following block- 
diagonal form on V <g> V with respect to the decomposition V = Vo © Va 

/ id \ 
id 
id 
V id J 

So, the block- diagonal form of H is permutation and, hence, obeys the relations 
( |4.4j) . Then the constructing blocks F^ act nontrivially on two sites i, j only, 
where they are the block-diagonal permutations ( f4.6|) . Now all the terms in ( |4.5| ) 
are pairwise interactions between two different sites and it can be written in 
explicit form 

dimV^ dimV^ 

E( V a v Q v v a\ , y y . \ ~* a \r a' 

\JL iQ A.j a + A ia A i0 J + A i0 Aj + 2^ A «a A ja' 



n 



HS 



E 



i 



" d 2 ■ 



a=l 



a,a '= 1 



The energy levels of 7i HS coincide with the levels of ( f4.1| ). The degeneracy 
degree with respect to the later is defined by (|3.14j) . The relations (|3.15| ), (|3.16| ) 
between the statistical sums of Hq S and H HS remain the same. 
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5 Fermionic representations 

In this section we consider the representations of some Hamiltonians derived in 
previous sections in terms of c ia , which are the annihilation-creation operators 
of spin cr =f, I fermion at site i. 

Below we are dealing with the quantum group U q sl2 only. We consider the 
simplest cases, where in the decomposition of V there is only one representation 
of spin > 0. Then i?-matrices, appearing in Hamiltonian are permutations and 
the Hamiltonian doesn't depend on q. So, it will be L^s^-invariant for all values 
of deformation parameter q. 

We identify the space V with the space of spin-a (er =t, I) fermionic wave- 
functions. Let us denote by |0) the fermionic vacuum: c a \0) = 0. Also we use 
|a> = c+|0>,UT> = -|U> = cJ-cf|0>. 

To construct the fermionic representations of Hamiltonians, introduced in the 
previous sections, it is convenient to use the fermionic representation of projection 
operators (|3.4|) , introduced earlier. Here we denote a, b = 0, ±1, 2, where means 
the empty site, ±1 mean the sites | |) and | j) correspondingly and 2 denotes 
the site | |f) with double occupation. Creation and annihilation operators can 
be expressed through the projection operators as follows 

c+ = X CT - aX 2 a C<7 = X° a -aX^ n a = X° + X\ a = ±1 (5.1) 

and vice versa 

X CT = (1 - n__ ff )c+ XI = wt. X°2 = c T q 

X a a = (1 - n_ CT )n CT X ° = (1 -71^(1 -n T ) = n h X\ = n^n i = d (5.2) 

Other formulae are obtained using (X^) + = X h a . Here we introduced the local 
Hubbard interaction operator d, which counts the double occupation of site and 
the hole number operator n h . 



5.1 t — J model for zero spin-spin coupling and Hubbard 
model in infinite repulsion limit 

Consider the simplest case, when the space of states decomposes into direct sum 
of fundamental spin-1/2 and trivial spin-0 multiplets V = V © V1/2 01 • We 
associate V with empty state |0) and V1/2 with | f), 
Using (|3.8|) for the present case, we obtain 



|) one electron states. 



L-l 



H(t,W x ,W 2 )=Y, 



i=l 



~t (^0^+1° + Xi+iyXi^) + H^iXjqX, 

a,a'=±l 





i+10 



+w 2 Yl x&x»ii 



a,a'=±l 
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After substitution of fermionic representation ( |5.2| ) of projection operators the 
Hamiltonian above (by canceling out nonessential in the thermodynamic limit 
boundary and constant terms) transforms into 



H t -j(t,W,n) = J E 

i=l 



-t 0- - U i -<r)(cfa C i+V + cf+l,aCi,a){l ~ n i+1 _ a ) 

a=±l 

L 

+Wnin i+ i} + fi^rii, (5.3) 

i=l 

where W = W\ + W2 and /i = —2W\. As TC(t,W,fi) conserves the particle 
number, it remains integrable without last constrain on /i. This Hamiltonian 
forbids double occupied sites and coincides with one of t — J model for vanishing 
spin-spin coupling J = 0. Integrability of this model and its correspondence to 
XX Z spin chain is well known |3"3fl . For vanishing density-density interaction 
V = (|5.3|) reduces to Hubbard model in infinite repulsion limit. 



5.2 Extended Hubbard with doping, pair hopping and 
hole density-density interaction 

Let consider now the case V = V © Vi. This corresponds to R\2{x) = R-2,i{x) = 
id in (|3.8|). Note that in this case H doesn't depend on Xi and q and we choose 
the coefficients there as follows 



H 



( Wi ■ id \ 

-t-id 

-t-id 

V W 2 ■ id / 



(5.4) 



We identify the empty state |0) with Vo and the space spanned by | |), | ||), 
I I) with Vi. The general Hamiltonian ( |3.9|) in our case has the form 



L-l 

H(t,W u W 2 ) = Yl 

1=1 



-t ^2 (.XiQX i+1 a + Xj + iQXj°) - t(Xilx i+1 



+Aj + ioAj2) + H^iXjgXj+XQ +W 2 ^ (Aj^ + Aj2)(Aj +1 ^, + A i+1 



<t,<t'=±1 

After substituting ( |5.2| ) into (|5.5| ) one obtain 



(5.5) 



L-l 



H(t,W 1 ,W 2 ) = Y, 



i=l 



t ■ - m _ <T )(cJ r c i+ i, (T + 4 +1>a c ita )(l - 

-*( c Jr c i c *+l.-l c «+i,T + c i+l,T c i+l,i c '-i c i,T) + ^l^NiVi ( 5 - 6 ) 

L 

+H/ 2 (n, - di)(n< + i - - 2H/ 2 £ n { + W 2 (n^ + n£) + W 2 L 



i=l 
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The quantum group generators ( |2.1| ) (after some renormalization) can be written 
in terms of fermionic operators also: 



ei 
fi 



[2]- 1 / 2 (n T c+ + n^) e = h 
[2]- 1/2 Kc|+^c+) f = ei 



(5.7) 



hi = 2 (nt 



n\ 



hr 



-hi 



The first term in Q5.6| ) is known as "doping" term, the second one is pair hopping 
term. Using n h = 1 — n + d and omitting unessential boundary and constant 
terms in (15. 61) one obtain 



L-l 



-* ZX 1 _ rii- a )(ct ifT c i+ i t<7 + c^ 1)CT c i)(J )(l - n i+ i t<7 ) 

(7 

t ( c i,1<$l c i+ 1 >l Ci + 1 >l + c 4 + +i,T c ti,l c *.i c i,T) + Wn i n i+i] (5-^ 



i=l 



where W = W x + W 2 , // = 2W 2 , U = -2W 2 . 

In view of the consideration, carried out in sec. 3, the energy levels of obtained 
fermionic model coincide with the levels of spin chain model 7io(t, W\, W 2 ) on W°. 
It is easy to see that its two site interaction term 



( W x \ 

-t 

-t 

V o o o w 2 J 



(5.9) 



leads to XX Z model in external magnetic field (see ( |5.15| ), (|5.16|) below). 

We note that the Hamiltonian ( |5.6|) preserves the number of holes Af h = 
J2i=i n i an d the number of double occupied states V = Y.f=i di- So, we can 
choose the parameters of corresponding terms in (|5.8|) to be arbitrary without 
loose of integrability. Of course, the Hamiltonian ( |5.8| ) is not invariant with 
respect to quantum group (|5.7|) now. 



5.3 Hubbard model with bond-charge interaction and ad- 
ditional density-density and boson-boson interactions 

Let us consider V(x a ,Xb) = V a (x a ) © Vb(x&) and identify the first multiplet in 
this decomposition (V a ) with the space spanned by | f), | I) and the second one 
(Vft) with the space spanned by |0), | ||). We consider two different commuting 

quantum group actions on V. The first one (U q sl 2 ^) acts on V a as spin-1/2 
representation (V a ~ V1/2) and on V& as two spin singlets (V5 ~ V © V ). The 

" ( 2 ) 

second one (U q sl 2 ) acts vice versa: V a ~ V © V and V b ~ V1/2. 
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Now by use of the form of intertwining operator H ( |3.5| ) for the decomposition 
V = V1/2 © V © Vq and taking vanishing matrixes A i and Ai Q , one obtain 



h(a,b) = y: 



EA S 1 ' V vS' | A V a 1 V <j3 



(7 % .<j3 



i/ 5 0|# A »<S< A i<7 + ^i05' A i CT A i# 



(5.10) 



Note that the indexes i,j here denote the nearest neighbors but not nonequivalent 
representations Aj, as in ( |3.5| ). So we used the upper index for them not to mix 
with ( p.5|) . Also, we changed the double indexes in (|3.5| ) on single indexes a — ±1, 
5 = 0,2by(l,ai)~cr, (no, 1) ~ S. Let us choose the parameters in ( |5.10| ) in the 
following way 



a <5*V 
n 00 S^Si 



8%5f 3 W 2 , Aii=W u B i 



5j 



R gi 



After these simplifications we have 



L-l 



H(t,W 1 ,W 2 ) = J2 



1=1 



~t (XigXi+il + X i+ x a s X ia ) 



(T,ct'=±1 
M'=0,2 



+ty a £ +w 2 x: 

ct,ct'=±1 <S,<5'=0,2 



(5.11) 



The nearest neighbor interaction of this Hamiltonian conserves its form with 
respect to index exchange 5 <-> cr. This means that 7"^(£, Wi, W2) is invariant with 

respect to both quantum groups U q sl 2 a,ndU q sl 2 ■ So, fl5.1l|) has U a sl 2 ®U a sl 2 - 
symmetry. The highest weight space at each site is formed by the two vectors 
I \) and |0) of V a and V 6 correspondingly. 

Using the fermionic representation of A-operators (|5.2p we obtain 



L-l 



H(t,W u W 2 ) = J2 



i=l 



+t(cf i0 .c i+li<7 + cf+^Ci^rn-a + n i+1 _ CT )} 
+ (Wt + W 2 )(n,jn i+1 - 2n,id i+l - 2n i+1 di + Adid i+ i) 
+ 2W l (d i + di+i) - W^m + n i+1 ) + Wi] 



(5.12) 



Two types of hopping are allowed here. First, fermion from single occupied 
side can hope to empty site and vice versa. Second, from double occupied site 
fermion can hope to single occupied site, which contains fermion with opposite 
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spin. Such hopping term appeared in [|TJ] after projection of Hubbard model 
on a fixed occupation number subspace. Thus, the Hamiltonian (|5.12|) preserves 
the occupation number of fermions: number of holes J\f h and number of double 
occupied sites T>. So, corresponding terms without affecting on the integrability 
an be rewritten with arbitrary coefficients. Also, note that the term 

L-1 

2(Wi + W 2 ) 2 did i+1 - riid i+1 - n i+1 di 
i=i 

can be written as hole and double occupation site interaction and Hubbard in- 
teraction in the following way 

L-1 

2{W l + W 2 ) £ n h t d l+1 + r$ +1 di -d { - d l+1 

i=l 

Thus one can generalize ( p,12| ) by 

L-1 

H{t,W u W 2 ,U,n)=Y. E [-*(<*c*+v + cki,a<v) 

i=l <r=U 

+t( C t<T C i+h<7 + C t+l,a C i,v)( n i,-v + ,-„)] (5.13) 

L-1 L L 

+W {ntm+x + 2n*d i+ i + 2nf +1 d 4 ) + U ^ di + V E "* 

i=l i=l i=l 

+W 1 (n 1 + n L ) + 2W 2 {d 1 + d L ) 

The term with coefficient W = W\ + W 2 consists of standard fermion density- 
density interaction between nearest neighbors and interaction between empty 
site (hole) and nearest neighbor site occupied by two fermions with opposite 
spins. The later interaction can be considered as interaction between two different 
bosonic sites |0), | {]). The boundary terms, which appeared in ( |5.13p , do not 
affect on the spectra of (|5.13|) in the thermodynamic limit and can be omitted. 
So, in this limit the Hamiltonian above depends only on W = W\ + W 2 . 
The connection between Hamiltonians ( |5.12| ) and ( [5. 13| ) is given by 

H(t,W 1 ,W 2 ,U,fi) = H(t,W 1 ,W 2 ) + (U+4W 2 )V+( f i+2W 1 )Af-W 1 (L-l) (5.14) 

Now we consider the restriction of ( |5.12| ) on the highest weight space, which 
is generated by the empty sites and sites occupied by the single a =\ spin. 
Looking on the matrix form of H (t, Wi, W 2 ) = A we recognize XX Z spin chain. 
It follows from our previous investigation that the obtained model is exactly 
solvable and has the same energy levels as XX Z Heisenberg model in the external 
homogeneous magnetic field B = Bz along z axes 

Ho(t,Wi,W 2 ) = H X xz(t,A,B) 

= -^E(^m + ^<i + A ^ + i + f^), (5-15) 
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A 



B = -(Wi - W 2 ) 



(5.16) 



where 

W 

~2t' t 

All energy levels have the same degeneracy. Recall that we identified the spin 
states |+), |— ) of ( |5.15| ) with the highest weight states | f) and |0) of V corre- 
spondingly. 

It is easy to check that the Hamiltonian ( |5.11| ) , ( |5.12 ) commutes with 77-pairing 



operators |14j, [15], [16], [17 

L 
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i=i 



i=i 



t=i 



A/"-L = ]T(A^-X, 



°1 



(5.17) 



i=i 



which generate s/2 algebra. Using this fact and ( |5.14|) one can compute the 
commutation rules of fl5.17| ) with generalized Hamiltonian (|5.13 ) 



[H(t, Wt, W 2 , U, ti),rj\ = -{U + 2^ + 4W)r) 
[H(t, W u W 2 , U, fi), V + ] = (U + 2» + AW) V + 



(5.1* 



So, we have [H(t, Wi, W 2 , U, //), r^ ± ] = if parameters satisfy the condition ji = 
-U/2 -2W. 

We are interested in thermodynamic behavior at T = 0. Denote by En = 
E N {t, Wi, W 2 , U, /z) the ground state energy of TC(t, Wi, W 2 , U, /i) in the iV-particle 
sector. Also, let En x = £jv x (A) is the ground state of XX Z spin chain ( |5.15| ) 
without external magnetic field (B = 0) in the sector with N± upturned spins. 
Recall that single occupied sites of particle space correspond to positive spin in 
XX Z model whereas double and empty sites correspond to negative spin. Thus, 
in the sector, which consists of states with fixed occupation number J\f\ — iVj., 
T> = (N — JVi)/2 the minimal energy Eft 1 of H(t, W u W 2 , U, pi) is (see (^Tg)) 



E 



iVi 
N 



£ 



u 



u 



— + W)N 1 +[ — + n + 2W)N+ -{W x - W : 



L 



(5.19) 



To obtain ground state En one should minimize (|5.19 ) among all possible values 
of JVi 



E 



N 



min Ep 



mm 

0<N!<N 



u 



L 



(5.20) 



+ [- + fi + 2W)N+ -(Wx - W, 



Note that the Hamiltonians TCxxz{t, A, 0) and Hxxz{—t, —A, 0) are related by 
similarity transformation |36[]. So, the spectrum of 7~Cxxz{t, A, 0) in the ferromag- 
netic regions (— t > 0, A < —1) and (— t < 0, A > 1) coincides. Hereafter we dial 
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with the first region and choose the parametrization A = — cosh 7, 7 > 0. Note 
that |U5D in zero magnetic field has two Z 2 symmetric ferromagnetic ground 



states |±, vac)xxz = |±, ±, • • • , ±) with energy £ = El = —t/2(L — 1)A. 

The energy levels of XXZ spin chain are determined by means of Bethe 



Ansatz p7| , 3SJ. In the thermodynamic limit L — > 00 there it has string type 
solutions. It was shown in [18| that the minimal energy in the sector with Ni 
upturned spins is 

sinhiVi7 „ , JVi , 

C7V1 = — 2tsinh7 — — + £ = — 2tsinh7tanh — 7 + £ (5.21) 

coshiv 1 7 + 1 2 

and corresponds to single string of length N±. Here the magnetization is restricted 
by < Ni/L < 1/2, i.e. we are dealing with excitations near ferromagnetic 
vacuum |— ,vac). Due to Z 2 symmetry we have £n 1 = £-l-Nx and the excitations 
near |+, vac) (1/2 < N\/L < 1) can be considered in a similar way by exchanging 
Nx -» L-N x in (^2ll). 



For infinite number of turned spins lini^oo Ni/ L = 0, L the minimal energy 



( 5.21 ) is Soo = — 2tsinli7 + 8q. So, all ground states for fixed magnetizations, 
which differ from both ferromagnetic ground states by infinite number of turned 
spins, have the same energy. 

It is well known that XXZ spin chain for |A| > 1 is massive. The energy 
gap between ground state and elementary excitations in the limit N\ — > 00 is 
AE = -2t(cosh7 - 1). 

Substituting ( p.21[ ) into ( |5.2U| ) one obtain the value of Ni, which minimize 



in (1O0D 



™ = /0 (U/2 + W<0) or (0 < U/2 + W< 7 /2, N < N c ) , 
1 \ N (U/2 + W>0) or (0 < U/2 + W < 7 /2, > N c ) , 1 > 

where N c is solution of 

N c fU \ U 
tanh -S 7 = - + W ) JV C , - - iy > 0, A^ c > 



2 ' V2 / ' 2 

Consider the regions (C//2 + < 0, < 2b) and (0 < U/2 + W < 7 /2, N < 
N c , W < 2t) in ( p.22|) , where there are no single occupied states. The ground 
state of Hamiltonian in the iV-particle sector (N is even) is any state with D = 
N/2 double occupied sites and L — N/2 empty sites. Let us write down all these 
states in the form 

^ = (ri+) N/2 \vac) , \vac) = |0, 0, . . . , 0) (5.23) 
Here r]p are generalized P-momentum //-pairing operators 

L L 



vp = Y, (i ' ,Jr j- r j-. vt = E e iPjc ti c h 

3=1 3=1 
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Note that t]q = r] ± . 

All the states ip£ obey ODLRO, i.e. 



lim ^ P '^X;^ P ^° ( 5 - 24 ) 



Note also that r/p operators obey rjp\vac) = 0. As a consequence the ground 
states ipp are superconducting. 

In the phase space regions (U/2 + W > 7/2, W < 2t) and (0 < U/2 + W < 
7/2, N > N c , W < 2t) the are no double occupied sites in the ground states of 
( |5.13| ). These states are linear combinations of states with N single fermions of 
same spin. Of course, they do not exhibit ODLRO Q5.24 ). 



The projection of ( |5.13| ) on states without double occupation (V = 0) gives 
rise to t — J Hamiltonian for zero spin-spin coupling J = ( |5.3D with additional 
boundary term W\(ni +til)- As it was mentioned above, this model is integrable 
and equivalent to Heisenberg magnet. Note that this equivalence can be seen 
from (|5.3|) . If we restrict T~C(t, W, fi) on the states with only one sort of fermion 
(spin- 1 or spin-j), we obtain the fermion representation of XX Z model 

L-1 L-1 L 

1<XXz(t, Wl, W 2 , /i) = _ *( C i" C i+l + C t+l C i) + W H U i n i+1 + ^Yl U i 

i=l i=l i=l 

For the special values of interaction potentials W\ = —W 2 = W the Hamilto- 
nian under consideration simplifies drastically. We have up to boundary terms 

L-1 

n(t,w) = J2 E [-*(<^ +v + c + 1)CT(V ) 

L L 

+t(cf i(7 c i+ i j(7 + cf +1)(T Ci t<r ) {ni- a + n i+ x- a ) +Uj2di + fi^rii (5.25) 

i=i i=i 

This is Hubbard model with bond-charge interaction. On the highest weight 
space it coincides with XY model in external magnetic field (or equivalently 
free fermion model with chemical potential term added) (A = in (|5.15j) ). So, 
TC(t, W) is exactly solvable and have the same energy levels as free fermion model 
with chemical potential term with degeneracy degree of each level. 

The Hamiltonian ( |5.25| ) T> was obtained and solved in [[12], [13]]. In our con- 



struction its integrability is obvious, because the restriction H°(t, W = 0) of 
( |5.25| ) on the highest weight space coincides with ordinary XY spin chain and, 



consequently, it is equivalent to free fermion spin chain. It was shown in |13[ that 
( |5.25| ) has superconducting ground states. Note that here we didn't consider the 
antiferromagnetic region, which corresponds to A = case. We suggest that 
in this region ( |5 . 1 3| ) has superconducting ground state too, obeying //-pairing 
mechanism. 
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